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Abstract 

We show that the conclusion on the breakdown of the standard small po- 
laron theory made recently by E.V. deMello and J. Ranninger (Phys. Rev. 
B55, 14872 (1997)) is a result of an incorrect interpretation of the electronic 
and vibronic energy levels of the two-site Holstein model. The small po- 
laron theory, when properly applied, agrees well with the numerical results of 
these authors. Also we show that their attempt to connect the properties of 
the calculated correlation functions with the features of the intersite electron 
hopping is unsuccessful. 

1 Jn a recent publication ]I| de Mello and Ranninger have analyzed numerically the familiar 
two-site Holstein model |2||| of a single electron coupled to an intra-site vibration mode. 
The model is the electronic doublet, <f>x, 2 , describing an electron legalized on sites 1 and 
2, respectively, plus the interaction with a vibration mode. The overlap of <f>\ and 02 leads 
to a splitting It of the doublet in the absence of the interaction, where 



t = J ftHafadV, (1) 
H e i is the electron Hamiltonian. The Hamiltonian of the model is given by 



p 2 MuJix 2 . + t t t 

= 2M H 2 _ a 2 a 2j - t{a[a 2 + a 2 ai). (2) 

Here p, x are the momentum and co-ordinate of the vibration mode, M, ujq its mass and 
frequency, respectively, and af, a$ (i = 1,2) are the electron operators. The Eq.(||) is an 
invariant with respect to the inversion x — > —x, (1, 2) — > (2, 1) and the parity is conserved. 
This model is an extreme simplification of the small polaron model. 2t is the analog of the 
electron bandwidth in a crystal. 

2. One of the basic results of the small polaron theory is the narrowing of the electron 

band due to the electron-phonon (vibron) coupling. In the strong-coupling regime, when 
g 2 —>■ oo, the narrowing is described by a simple exponential dependence on g as 

t* = texp(— const • g 2 ). (3) 

The exponential dependence, Eq. (2) is readily derived by the use of the familiar double- well 
potential model f2|. Polaron tunnels coherently within the narrow band at low temperatures 
(while the thermally activated hopping dominates at high temperatures HH). The numerical 
solution for several vibrating molecules coupled with one or two electrons |§Hl2| revealed 



an excellent agreement of the numerical bandwidth with the analytical Holstein and Lang- 
Firsov (LF) results at large g. 

3. However, de Mello and Ranninger arrived with an opposite conclusion. On the basis 
of numerical analysis of the same problem authors [Q claim in Sec. Ill that: 

the LF approach^, which is generally believed to become exact in the limit of 
antiadiabaticity and an electron-phonon coupling going to infinity, actually di- 
verges most from the exact results precisely in this limit... (p. 14885) 



1 In Ref. (T| the "LF approach" is identified with the lowest (~ t) order in perturbation theory. 
We do not agree with this identification, since to go beyond the lowest approximation was a central 
point for the authors of Ref. [UH0] . 
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These authors did not provide any physical explanation for their drastic disagreement 
with all earlier results starting from the pioneering work by Holstein and including the 
kinetic theory of strongly-coupled electron-phonon systems, in particular with the theory of 
high-frequency conductivity [0,0]. 

4. In Ref. |TJ[ we have recently developed the analytical approach to the two-site model by 



the use of the expansion technique, which provides the electronic and vibronic terms as well 
as the wave functions in any order in powers of t. In the second order in t the doublet energy 
of the ground state, E± is given by 

E ± = ±f-^ ) t* = texp(-2E p /huj ). (4) 

Here E p = g 2 /2Mu 2 is the polaron shift (the following designation is used in Ref. [0: 
E p = a 2 huo). The first term describes the splitting of the doublet (components of the 
doublet have opposite parity) corresponding to the bandwidth in a crystal, as discussed 
above, while the second term is a correction to the polaron shift of the whole band due to 
the virtual transitions to the nearest-neighbour site. The exponential reduction factor was 
found in all orders of t of the perturbation expansion [|nj] in agreement with the standard 
result, Eq.(2). On the other hand the corrections to the atomic level are relatively small as 
l/g 2 rather than exponential. 



5. In Ref. ||14j| , Sec. 6 it was demonstrated, that the aforesaid statements of authors of 
Ref. Jl| do not correspond to reality and is only due to the fact that they failed to notice 
the above-mentioned difference between the splitting of the doublet's components and their 
shift as a whole. This is the result of the methodological defect of the approach which was 
employed in Ref. |IJ. Instead of a direct solution of the quantum-mechanical problem (to 
determine the energy spectrum and the wave functions []), they calculated a value (which 
has no direct physical meaning) 

£ ki n ~ t{a\a 2 ) h (5) 



2 This approach would eliminate the possibility to make such an error to a considerable extent. 
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(. . (i = +, — ) is a quantum-mechanical average on the one of the doublet's component. 
They accepted implicitly an assertion that: 

El in is the analog of the electron bandwidth 



As it was shown in Ref. [H|, this assertion is incorrect. Let us denote as AE l a correction 
term to the energy level % which is generated by the last term ~ t of the Hamiltonian (2). 
In fact: 

1. El in is proportional to dAE l /dt (rather then to AE l ). The average (5) decreases when 
\g\ increases as a power of g 2 (~ g~ 2 for g 2 — > oo, see Ref. Il5|), but not exponentially! 



2. Analog of the electron bandwidth is the difference SE = \ AE + — AE |, but not AE l 
alone. SE contains the exponential factor (3). 

In other words, in Ref. authors had compared the values of the essential different nature. 
This is the source of the above-mentioned drastic disagreement. 

6. In Ref. [I]], Sec.V authors have calculated the electron dipole momentum^ and the vibronic 
co-ordinate correlators (deformation dynamics correlator in Ref. [|1|]): 

Xnn(r) = ((m - n 2 ) T (n 1 - n 2 )) , Xxx(r) = (x T x) , (6) 

where (. . .)o is the average over the ground state \l/o- Not any analytical examination was 
made. The authors of Ref. Il|] pointed out that the calculated curves which represent a 



3 No explicit wording of this assertion is given in Ref. However, this wrong interpretation of 
££ in is forced, because the authors of Ref. |TJ] have compared it just with a small polaron bandwidth 
I^Q. In other case such comparsion would be irrelevant, since authors of Ref. HQ never calculated 

-^kin- 

4 The authors of Ref. (l[ have connected this correlator with the charge fluctuations which is not 
quite correct. 
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functional dependence x( T ) (Figs- H> 12 of Ref. |IJ) may be presented as a superposition of 
a slow and fast oscillations with the frequencies of the fast oscillations t and u> accordingly 
for Xnn(T) and Xxx(t) Q- They claim that when these frequencies (which are certainly > Uo, 
see Table 1 in Ref. [[![]) draw together, the qualitative changing of the electron transport 
mechanism takes place. No physical argumentation (even of a qualitative nature) to support 
this assertion of the authors have been given. 
7. We note that x( T ) ma y De represented as 

X(r) = J2 a 2 m e- iuJm0T , hu m o = E m - E , a m = (^ m A^ ), (7) 

ty m , E m are eigen- functions and eigenvalues, accordingly, of the Hamiltonian Eq.@, A = 
rii — n 2 or x . The summation In (|7|) is performed over the states \l/ m with parity opposite 
to the parity of the ground state \l/o (selection rules for the operator A). We note that the 
frequency spectra of Xnn and Xxx are identical. 

There is a connection between correlators (6), (7), and corresponding generalized sus- 
ceptibilities K a (u) (see Ref. Here uo is the frequency of an external disturbance. For 
example, the complex polarizability k(uj) of the considered model may be expressed by the 
Fourier-transform of the correlator (6) 

k (u;) = —J q e^ + ^ (xna(r) - xU~t)) dr, 5 > 0, 5 - 0. (8) 

/ is a constant with dimensionality of the length. An imaginary part of k(u) is 

ie 2 l 2 

k ' (u)) = ——^a 2 m (5(uj - u m0 ) -S(u + u m0 )). (9) 

The value lok ' (uj) determines an absorption coefficient of electromagnetic radiation, and 
the value io m a^ m determines an absorption intensity for the transition — > m. Due to the 
presence of the 5— functions in Eq.(j^), the absorption process in the given frequency range 

5 We do not see any reason to identify t with " renomalized intrisinc hopping integral t" . 
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can not be linked causally with another one in the other frequency range []. 

It is natural to identify the frequencies in Ref. [IJ, Figs. 11, 12 in the following way: the 
slow oscillation corresponds to the width 5E = It* of the lowest doublet, Eq. (|4]); the fast 
oscillations correspond the frequencies > SE/h in Xnn(T) and Xxx(t) for that the weights o? m 
(see Eq.([7p) are maximal. The slow and fast oscillations are located in different frequency 
regions, therefore changes in the high-frequency region cannot modify the low-frequency 
electron transport mechanism cardinally. 

For these reasons, the aforesaid assertion (see above, Sec. 6) is unfounded. And again, 
the source of this error is of methodological character. This situation would be excluded, if 
instead of x( r )> which have no direct physical meaning, they have considered the complex 
polarizability k(u>). 

In Ref. Jl| the authors touch upon the subject of a boundary where a localized regime 
changes to an itinerant one. In our opinion the boundary is determined by the parameter 
rji = t/2E p []. For rji < 1 the lowest adiabatic potential curve has two minima, which 
are separated by the energetic barrier (localized regime), for rji > 1 the barrier vanishes 
(itinerant regime), Ref. H]. We state, that in the range of the parameters, which considered 
in Ref. |TJ the itinerant regime was not yet realized. 

8. We have checked and proved that (under right interpretation, naturally) the numerical 
calculations presented in Sec. Ill and V of Ref. fTI] agree satisfactorily with the Holstein-LF 



6 For example, a DC conductivity of the semiconductor does not depend on a higher empty band 
contribution (and vice versa, the intensity of the interband absorbtion does not depend on DC 
conductivity mechanism) . 

7 r]i is the parameter which was introduced by T.Holstein; the small polaron appears when r}\ < 1. 
Also this parameter determines a correction (~ t) to the overlap integral between the site-localized 
functions. 
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approach for t/hu < 10. 

Finally, we note that the authors' assertion in Sec. V [|TJ 

'We notice that the charge dynamics qualitatively tracks globally the behavior 
expected on the basis of the LF approximation in the anti-adiabatic limit... 
'(p.14882) 

obviously clashes with their statement in Sec. Ill (see above, the quotation in our sec. 3 ). 
No comments on this discrepancy are given in Ref. [[J . 

In conclusion we state that although the numerical calculations in Ref. [|lj were performed 
fairly enough, yet their interpretation is untenable. 

One of us (V.V.K.) acknowledges support of the work by RFBR Grant 97-2-16705 and 
the Ministry of Science and Technology of Slovenia. 



8 The deviations which arise when t/hu)Q > 1 may be explained qualitaively in the framework of 
the adiabatic approach. 
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